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ABSTRACT 

We observed the Crab pulsar with the 43-m telescope in Green Bank, WV 
over a timespan of 15 months. In total we obtained 100 hours of data at 1.2 GHz 
and seven hours at 330 MHz, resulting in a sample of about 95000 giant pulses 
(GPs). This is the largest sample, to date, of CPs from the Crab pulsar taken 
with the same telescope and backend and analyzed as one data set. We calculated 
power-law fits to amplitude distributions for main pulse (MP) and interpulse (IP) 
Q_f CPs, resulting in indices in the range of 2.1—3.1 for MP GPs at 1.2 GHz and 

in the range of 2.5-3.0 and 2.4-3.1 for MP and IP GPs at 330 MHz. We also 
correlated the GPs at 1.2 GHz with GPs from the Robert C. Byrd Green Bank 
Telescope (GBT), which were obtained simultaneously at a higher frequency (8.9 
GHz) over a span of 26 hours. In total, 7933 GPs from the 43-m telescope at 1.2 
GHz and 39900 GPs from the GBT were recorded during these contemporaneous 
observations. At 1.2 GHz, 236 (3%) MP GPs and 23 (5%) IP GPs were detected 
at 8.9 GHz, both with zero chance probability. Another 15 (4%) low-frequency 
IP GPs were detected within one spin period of high-frequency IP GPs, with 
a chance probability of 9%. This indicates that the emission processes at high 
and low radio frequencies are related, despite significant pulse profile shape dif- 
ferences. The 43-m GPs were also correlated with Fermi 7-ray photons to see if 
increased pair production in the magnetosphere is the mechanism responsible for 
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GP emission. A total of 92022 GPs and 393 7-ray photons were used in this cor- 
relation analysis. No significant correlations were found between GPs and 7-ray 
photons. This indicates that increased pair production in the magnetosphere is 
likely not the dominant cause of GPs. Possible methods of GP production may 
be increased coherence of synchrotron emission or changes in beaming direction. 

Subject headings: Crab pulsar, Fermi, giant pulses 



Introduction 



The Crab pulsar was discovered by Staelin & Reifenstein in 1968 through its giant 
pulses. Giant pulses (GPs) can be thousands of times brighter than the average pulse. The 
temporal occurence of GPs is random but at frequencies below 3 GHz they always occu r at 
the phase of either the Crab pulsar's main pulse (MP) or interpulse (IP) (lLundgrenlll994[ ) . At 
frequencies from 4 GHz to 8.4 GHz, GPs are emitted at the phase s of the MP and IP as we ll 
as at the phases of two additional high-frequency components (IMoffett fc HankinsI 119961). 
Abov e 8.5 GHz, GPs are again seen at only the phases of the MP and IP (jJessner et al. 



2OIOI) 



The emission mech anism of GPs is still an open question (e.g. iPetroval l2006l : llstomin 
2OO4J : IWeatheralll Il998l ). GPs could be caused by increased pair production in the mag- 



netosphere, increased coherence of synchrotron emission, or changes in beaming direction. 
Correlating radio GPs with high-energy photons is one way to deterra ine if increased pair 
production is a major cause of GPs. A recent model (ILyutikovl 120071 ) proposes that GPs 
are generated on the last closed magnetic field line near the light cylinder via anomalous 
cyclotron resonance. If this is true, there would be an increase in 7-rays at the times of 
radio GPs. An omalous cyclotron res onan ces have been previously proposed as a cause of 
GP emission by iLyutikov et al.l (Il999[ ) and iMachabeli fc Usovl (119791 ). We expect the 7-rays 
to be phase aligned with the radio GPs at 1.2 GHz as the 7-ray and 1.2 GHz radio profiles 
are aligned (see Figured]). Since we have such a large data set of GPs, we can correlate 
them with 7-ray photons to test this mod el. Correlations be twee n radio GPs and 7-ra y 
photons have been c arried out previously by lBilous et al.l (1201 ll ) and lLundgren et al.l (119951 ). 
Shearer et al.l (l2003l) correlated r adio GPs and optical photons and found a slight correlation, 
and work by ICoUins et al.l (|2012| ) supports this result. 



The Fermi Large Area Telescope (LAT) is a pair conversion telescope that operates in 
the energy range from 20 MeV to 300 GeV. It has a large field of view (2.4 sr), which allows 
it to rapidly map the entire sky, and very good angular resolution, minimizing background 
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contamination (lAtwood et al.ll2009r). The Cra b pulsar is a bright Fermi source and has been 
studied at these energies. lAbdo et ahl (|2010[ ) found that the 7-ray profile is double peaked 
and matches the 1.4 GHz profile, with the 7-ray peaks leading t h e rad io peaks by ~0.01 
phase. Since the LAT covers the energies predicted by iLyutikovl (120071 ) and observes the 
Crab pulsar multiple times per day, we can use Fermi data to test Lyutikov's theory. 



The power-law nature of t he amplitude distribution of CPs is well known (e.g. lArgyle fc Cower 



I972I : iPopov fc StappersI 120071 ) , but there hav e been varying values for the power-law inde x 



calculated, even for similar frequencies (e.g. iKaruppusamy et al.l I2OIOI : iBhat et al.l l2008l ). 
With our large sample of CPs, we can calculate power-law indices over a long tirnespan . 
Since the amplitude distribution of pulses from many pulsars is log- normal ( lRitchingslll976l ) . 
CPs must have a different emission mechanism. By comparing our power-law index with 
those of other neutron star source classes (pulsars, magnetars, RRATs), we may be able to 
put constraints on the CP emission mechanism. 

A 7-ray flare was recently observed from the Crab Nebula by the ACILE satellite 



(iTavani et al. 



2011 ) . Th e flare lasted from MJDs 55457—55461 and was also observed by 



Fermi ( iHays et al.l I2OIOI ). This occured during the span of our observations, so we can 
compare radio and 7-ray properties of the pulsar from before the flare with a few days of 
data taken about two months after the flare to verify that the flare was not associated with 
the pulsar. Current theories suggest that the flare was caused by electro ns accelerated by 



magnetic reconnection at the termination s hock of the Crab Nebula (e.g. iBednarek fc Idee 



2011; 


Bvkov et al. 


2011; 


Cerutti et al. 


2011) 



The plan for the rest of this paper is as follows. Section[2]describes the radio observations 
and outlines the radio analysis. Section |3] discusses the CP statistics and the power-law 
index calculations. Section H] details the correlation between 43-m and CBT CPs. Section 
E] presents the Fermi data used and its reduction as well as the 43-m/Fermz correlation. 
Section [6] discusses the recent 7-ray flare. Finally, conclusions are offered in Section [3 



2. Radio Observations 



We collected data over the course of 15 months using the 43-m telescope (hereafter 
CB43), located at the National Radio Astronomy Observatory site in Creen Bank, WV. 
The CB43 is funded by M IT Lincoln Labs to perform bistatic radar observations of the 
ionosphere ( ILangstoru 120071 ) . Most of the observations had a center frequency of 1.2 CHz 
with a useable bandwidth of 400 MHz, with a handful of the later observations centered at 
330 MHz, with a useable bandwidth of 150 MHz. All of the observations were taken using 
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4096 frequency channels and varying sampling times, listed in Table [T] Some of the sampling 
times used were rather long, and this was done to decrease the data volume. In total, we 
observed the Crab pulsar for 100 hours at 1.2 GHz and seven hours at 330 MHz. The 330 
MHz mode was implemented later and used b ecause pulsars, in gen eral, are stronger at lower 
frequencies and exhibit a steep spectral index ( iLorimer et al.lll995l ). Even though the smaller 
bandwidth at 330 MHz reduce s sensitivity by alni ost a factor of 2, the steep spectrum of 
the Crab pulsar (a = 3.1±0.2 (ILorimer et al.l 119951 ). where S = means that the pulsar 
is 55 times brighter at 330 MH z than at 1.2 GHz. The flux density from the Crab Nebula 
is greate r at lower frequencies (ICordes et al.l |2004| ) but, due to a shallow power-law index 



-0.27\ 



flAllenlll973l : iBietenholz et al.lll997n . at 330 MHz (flux density ~1288 Jy) is only 1.4 



times brighter than at 1.2 GHz (flux density ~909 Jy). Therefore there should be many 
more CPs above our lOo" threshold at 330 MHz than at 1.2 GHz. 

Data were taken with WUPPI, the West Virginia University Ultimate Pulsar Processing 
Instrur nent. WUPPI is a c lone of GUPPI (Green Bank Ultimate Pulsar Processing Instru- 
ment) (IRansom et al.ll2009l ) for use with the GB43. GUPPI is a flexible digital backend for 
the GBT. Like GUPPI, WUPPI is built from reconfigurable off-the-shelf hardware and soft- 
ware available fr om the CASPER (C enter for Astronomy Signal Processing and Electronics 
Research) group (jParsons et al.ll2009l ). Both GUPPI and WUPPI sample the data with 8-bit 
precision over bandwidths as large as 800 MHz, and are capable of recording all four Stokes 
parameters. To ease disk space usage, only total intensity (Stokes I) was recorded. 

The data were processed using a real-time data reduction pipeline on a 16 processor 
mini-cluster. The pipeline is comprised of a set of scripts, built from freely avail able analysis 
softwar^, that can reduce a file in real time. First, a file in PSRFITS format ( Hotan et al. 



20041 ) has its mean bandpass divided out. Frequency channels that have intensities above the 
resulting mean are flagged as containing radio frequency interference (RFI) and removed. 
Then the resulting file is dedispersed at both a dispersion measure (DM) of zero and the DM 
of the Crab, which is ~56.8 pc cm~^; the exact val ues for each obse rvation, obtained from 
the Jodrell Bank Crab Pulsar Monthly Ephemeri^ ( Lyne et al. 1993), are listed in Tabled] 
The DM of the Crab pulsar can vary on a monthly timescale by ~0.01 pc cm~^. In order 
to keep processing in real-time we did not correct for this. To ease space requirements, raw 
data were not kept, so we could also not correct for this later. Errors in CP arrival times 
due to an incorrect DM on the order of 0.01 pc cm^^ are ~20 fis, which is less than the 
sampling time. Both resulting time series are then searched for single pulses with a peak 



http : //sigproc . sourcef orge .netl 



' http : //www . j b . mail . ac . uk/ -pulsar/crab . html 
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signal-to-n oise (S/N) above lOcr using our single pulse search (based on the algorithm de- 
scribed by lCordes &: McLaughlin! (j2003[ )). Our 10a definition of a GP is somewhat arbitrary 
because there is no set threshold for a GP as the weakest GPs have yet to be observed. The 
current population of observed GPs accounts for no less than 50% of the pulsed emission at 
frequencies around 1.6 GHz, and the inclusion of the weakest GPs could bring that number 
up to 90% (IMajid et al.ll201ll ). A GP time-of-arrival (TOA) (measured by taking the arrival 
time of the peak of the pulse), the S/N of that peak, and the pulse width, taken as the width 
of a best-fit top hat function, are recorded for each pulse. The TOAs from each time series 
are compared and if a pulse detected at the DM of the Crab is within 0.1 s of a zero DM 
pulse, the S/N is checked for both pulse^. If the S/N is higher in the zero DM pulse, the 
pulse is assumed to be due to RFI and removed. The pipeline outputs a profile for each 
GP and produces an average folded profile for the entire observation, a plot of GP arrival 
time vs pulse phase, and an average GP profile, made by summing all of the individual GP 
profiles. 

Figure [2] shows the number of GPs detected versus pulse phase. It is worth noting that, 
although there were only four post-fiare observing epochs, there are many more GPs than 
were recorded in all of the pre-fiare data. This was due to a receiver upgrade which resulted 
in a large increase in sensitiv ity between our ob servations on MJD 55412 and MJD 55516. 
At frequencies below 5 GHz ( ICordes et al.l 120041 ) . most detected GPs come at the phase of 
the MP, with the rest at the IP phase. In our data, ~87% of detected GPs were at the phase 
of the MP, while th e other ~13% were at the phase of the IP. Only 5% of GPs detected by 
Cordes et al.l (120041 ) at 1.2 GHz were at the phase of the IP, but iKaruppusamy et al.l (l2010l ) 
found that ~12% of GPs detected at 1.4 GHz were IP GPs. For detected pulses that show 
up out of phase with the MP and IP, the frequency versus time plots are checked by eye to 
see if they show the proper quadratic frequency sweep for a DM of 56.8 pc cm~^. If they 
do not, those pulses are removed. In most cases, the RFI removal process mentioned above 
removes most of these false detections before they are checked by eye. We found no events 
at phases other than those of the MP or IP that showed the proper frequency sweep to be 
a real GP. 



Most of the GPs we detected have a constant intensity across the entire band. Some 
of them, however , show variations in amplit ude as a function of frequency, which have been 



seen before (e.g. IKaruppusamy et al.ll2010l ). We also see variations in amphtude between 



days, which is likely due to refractive interstellar scintillation (RISS), which can affect the 
strength of pulses on a timescale of days. 



■^This is of the order of the frequency dependent arrival time delay across the band, and is used as it 
removes very nearly all RFI, while removing few real GPs. 
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In order to determine if the time stamps for the GB43 are accurate, a short (~10 
min) observation of the Crab pulsar was made contemporaneously on MJD 55406 at similar 
frequencies with the GB43 and GBT. The GB43 used a center frequency of 330 MHz, while 
the GBT used 350 MHz. GPs were detected in each data set and their TOAs were compared. 
As can be seen in Figure [3l the pulse shapes of contemporaneous GPs are similar, and the 
TOAs, converted to infinite frequency at the solar system barycenter, have identical arrival 
times to within the instrumental resolution. 

Figure H] shows the average profile of the Crab pulsar at 1.2 GHz and 330 MHz. The 
intensities are in arbitrary units. The MP and IP are visible at both frequencies. The 
increased width at low-frequency is due largely to DM smearing, which is on the order 
of 2.4 ms at the bottom of the band. The 1.2 GHz profile exhibits a weak low-frequency 



component (LFC), about one-tenth of a pulse phase ahead of th e MP fe.g.lMo: 



19961 ). At 330 MHz there is a precursor (PR) to the MP (e.g. iRankin et al. 



ett &: Hankins 



1970h which is 



broa d enough at this frequency that it shows up as the first peak of the double-peaked MP 



see 



Karuppusamy et al.l 120121 . for the frequency evolution of the PR), which is visible in 



the low- frequency GBT profile from MJD 55406 (Figured]). Due to a much larger effective 
resolution, the MP in the GB43 profile is unresolved and appears to have a single peak. 
Table [2] lists the effective resolutions at the top and bottom of the band for both the 330 
MHz and 1.2 GHz GB43 observations as well as the GBT 350 MHz observation. 

The full-widths-at-half-maximum (FWHMs) of the MP and IP from the GB43 profiles 
are ~305 /xs and ~360 fis at 1.2 GHz and -3.2 ms and ~1.5 ms at 330 MHz. At 330 MHz, 



the MP is wider than the IP due to its overlap with the PR (IRankin et al.lll970l ). The average 
GP profiles, however, are narrower than the average folded profiles at both frequencies (see 
Figure [5D, with FWHMs of ~140 ^s for both the MP an d IP at 1.2 GHz, an d ~1.1 ms at 330 
MHz for both the MP and IP. It has been suggested by iPopov et al.l (120061 ) that every pulse 
at the phase of the MP and IP is a GP, and that normal emission only comes from the PR, 
where no GPs have been seen. We also have not seen any GPs at the phase of the PR at 
330 MHz in our data. Even though our observe d GP profiles are narrower than the average 
folded profiles, this hypothesis is not ruled out. iPopov fc Stapperd (120071 ) found that pulse 
width is inversely proportional to pulse intensity, so lowering our S/N definition of a GP 
would include wider GPs which could possibly increase the width of the average GP profile 
to that of the folded profile. Also, the phases of weaker GPs may have more deviation from 
the center phase of the MP and IP than stronger GPs, leading to a wider average profile. 
Unfortunately, we were unable to lower our GP threshold below lOcr as this resulted in many 
spurious pulses. 
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3. Amplitude Distributions 



A total of 93698 GPs were observed with the GB43. At 1.2 GHz we recorded 76707 
GPs at the phase of the MP and 10871 GPs at the phase of the IP, and at 330 MHz 5232 
MP and 888 IP GPs were recorded. Due to low GP statistics on some days, only 77122 MP 
GPs and 2135 IP GPs were used in fitting the amplitude distributions. 

We calculated, through least-squares fitting, power-law indices for the differential ampli- 
tude distributions of MP and IP GPs at both 1.2 GHz and 330 MHz for each day separately. 
Figure |6] is an example of a power-law fit for the MP. We noted t h at the power-law index of 
GPs varied daily, as was previously reported by iLundgren et al.l (119951 ) at similar frequen- 
cies. This gave a range of power-law indices from 2.1±0.3 to 3.1±0.2 for MP GPs at 1.2 
GHz, 2.4±0.4 to 2.81±0.03 for IP GPs at 1.2 GHz, 2.5±0.2 to 2.95±0.09 for MP GPs at 330 
MHz, and 2.4±0.2 to 3.1±0.2 for IP GPs at 330 MHz. Table |3]hsts the power-law indices 
for each day. These indices agree with other published power-law indices, listed in Table 
H The rang e of in dices is likely due to wider GPs having steeper spectra, as was seen by 



Popov et al.l ( I2OO8I ). We did not have the time resolution to separate the GPs by width, as 



the widest GPs reported by iPopov et al.l ( l2008l ) were 64 /_fs, which is the shortest sampling 
time we used. 

We were unable to properly calibrate the data due to lack of any off-source pointings 
and/or observations with a pulsed cal. There were also doubts as to the stability of the 
receiver over long timespans. We were able to do a rough calibration, however, by correcting 
for RISS, which affects the str ength of all GPs by th e same amount, so there should be no 
change to the power-law slope. iRickett fc Lynd (Il990[ ) found that the RISS timescale for the 
Crab pulsar scales as Based on this, the timescales at 1.2 GHz and 330 MHz are 1.7 



days and 30 days, respectively. A Lomb-Scargle analysis (IScarglelll982l ) of the periodicity of 
GP arrival times yielded periodicities of 0.41 days at 330 MHz and 0.99 days at 1.2 GHz. The 
~one day periodicity of the GP arrival times in the 1.2 GHz data is on the order of the RISS 
timescale at that frequency, but the same periodicities were also seen in the randomized data 
set, so they are not significant. We attempted to correct for RISS by scaling observations 
based on the brightness of t heir folded p r ofiles . However, if GPs dominate the average 
folded profile, as was seen by iPopov et al.l (|2006[ ). any intrinsic variations in their strength 
will affect this correction. To determine if this was true, we removed the GPs from the 
average profile. We found that at 1.2 GHz the amplitude of the average profile without 
GPs was 1% less than the average profile with the GPs, and at 330 MHz was 2% less. We 
therefore concluded that GPs do not dominate the average p r ofile. One possible reason why 
our results are the opposite of what was seen by lPopov et al.l (120061 ) is that they recorded all 
GPs above 5a, while we were only able to record GPs above 10a. Therefore, it is possible 
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that weaker GPs dominate the average profile. If the average profile is truly dominated by 
GPs, then correcting for RISS is complicated, but since we do not see this, we used the 
following correction. 

We took the brightest folded profile from each of our three observing epochs (pre-flare 1.2 
GHz, 330 MHz, post-flare 1.2 GHz) and scaled all of the folded profiles from those epochs to 
that profile. That provided us with a scaling factor for each day, by which we then multiplied 
the GP S/Ns. We also compared our indices with power- law indices of other source classes 
in an attempt to constrain the GP emission physics. Other source classes and their indices 
are listed in Table [51 Although most normal pulsars have log-normal distributions, some 
pulsars have amplitude distributions that have power-law tails. These power-law exponents 
are included in Table under 'Normal Pulsars'. Our power-law indices match those of 
magnetars and RRATs, but do not match the power-law tails seen from normal pulsars. 

Also seen in Figure [6] is a non power-law tail, which is see n on all days in bo th the 



MP and IP. This deviation is significant, and was also seen by ICordes et al.l (120041 ). who 
postulated these outlying GPs could be supergiant pulses, indicating that there may be two 
distinct mechanisms for GP generation. These supergiant pulses account for slightly less 
than one percent of the GPs used in the amplitude distributions. 



4. Correlating Radio GPs from the GB43 and GBT 

For 16 hours over the span of eight days, we observed the Crab pulsar simultaneously 
with the GB43 and GBT. The center frequency for GB43 observations was 1.2 GHz with 400 
MHz of useable bandwidth, while that of the GBT was 8.9 GHz with 800 MHz of useable 
bandwidth. At high frequencies the majority of GPs come at the phase of the IP, with 
fewer MP GPs and few GPs from high-frequency components (HFCs) (see Figure [7] for a 
comparison of the folded profile at high and low frequencies). For this work, we did not use 
any HFC GPs in our correlation analysis. We matched the barycentered arrival times of 
the 39900 GPs with S/N > 10 recorded with the GBT (1035 MP GPs, 38865 IP GPs) with 
the 7933 GPs recorded simultaneously with the GB43 (7466 MP GPs, 467 IP GPs). We 
found that 236 low-frequency MP GPs were also detected simultaneously at 8.9 GHz. These 
GPs were neither the strongest nor the weakest pulses from either data set. The chance 
probability of this occuring is zero percent. This is not surprising, as we know that the high- 
and low-frequency MP are the same component. We did expect, however, that all of the MP 
GPs detected at 8.9 GHz would be the strongest of the MP GPs detected at 1.2 GHz, since 
the MP is significantly weaker at 8.9 GHz than at 1.2 GHz. All chance probabilities were 
calculated assuming Poisson statistics, using the formula P = where A is the number 
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of detections expected and K is the number of detections actually recorded. 



Moffett fc Hankind (119961 ) first noticed that the low- and high-frequency IP components 
of the folded profile are separated by 10°, which is about 970 /is. We checked to see if 
any low-frequency IP GPs were simultaneously detected at 8.9 GHz. Only 23 were found 
(~5%), with a z ero percent chance occurence, but this is expected since the spectral index 
of GPs is steep (jPopov et al.ll2008l ). We then checked to see if low- and high-frequency IP 
GPs cor nmonly occu r ed wi thin the same rotation period of the pulsar, as was previously 
seen by iPopov et al.l (120081 ). In our data, we found 15 instances (~4%) of low- frequency 
IP GPs occuring within one spin period of a high-frequency IP GP. The probability of this 
measurement happening by chance is 9%. In four of these instances the GPs are within 
one /US of each other and the chance probability of this occuring is zero. In the other 11 
instances, the GPs are almost one spin period apart. These are likely statistical, as there 
is a 10% chance of this occuring randomly. The fact that only four of the 1.2 GHz IP GPs 
occur within one fis of an 8.9 GHz IP GP s uggests that t he high and low-frequency IP 



may be created by different physical pro cesses (IMoffett 



19971') p ossibly due to emission from 



lyvypc 

different regions in the magnetosphere (IHankins fc Eilekl 120071 ). However, since the chance 
probability of detecting four low- and high-frequency IP GPs within one fj,s is z ero, it seems 



l ikely that the high-frequency IP is related to the low-frequency IP, as seen by iPopov et al. 



( 120081 ) in previous correlations of IP GPs at 600 and 4850 MHz. One possibility is t hat both 
the h igh and low-frequency IP are reflections of the MP off of the magnetosphere (jPetrova 
20091 ). In this scenario, radio emission from the MP originating deep in the magnetosphere 
propagates through the electron— positron plasma that fills the magnetosphere. Tranverse 
scattering causes the MP emission to be backscattered, causing it to arrive at a different 
pulse phase. If this was the case, we would expect a linear relationship between the strength 
of the MP and IP in the folded profile for each observation. As shown in Figure El we do 
see such a relationship at both 1.2 GHz and 330 MHz. The slopes for the 1.2 GHz and 330 
MHz relations are 0.310±0.004 and 0.58±0.08, respectively. 



5. Correlating Radio GPs from the GB43 with 7-ray Photons from Fermi 



One of the predictions of iLyutikovl (120071 ) is that there would be increased 7-ray flux 
during a GP. In order for this model toaccurately reproduce the data, the plasma density 
of the GP emission region must be ~10^ higher than the minimum Goldreich- Julian density 
and the duty cycle of the pulsar must be 0.001. This increase in density could be due 
to enhanced pair production in the pulsar magnetosphere. If correct, these high-energy 
particles, produced during reconnection close to the Y point, where the last closed magnetic 
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field lines approach the light cylinder at the magnetic equator (lLyutikovll2007l ). are expected 
to produce curvature radiation from 0.1—100 GeV (depending on the value of the Lorentz 
factor 7) at the time of a GP. This curvature radiation would cause an increase in 7-rays at 
the times of radio CPs. 

Data were downloaded from the Fermi online archiv^ for days when radio observations 
occured. 'Source' class events from Pass 7 data above 100 MeV were selected in an energy 
dependent radius {9 < Max(6.68 — l.TGlo gfElJ.S)", where E is the energy of the photon 



in MeV) around the position of the pulsar ( lAbdo et al.ll2010l ). Only photons in Good Time 



Intervals (GTIs) were selected and those with a zenith angle > 100° were excluded to dis- 
criminate against 7-rays generated in the Earth's atmosphere. The photons were converted 
to infinite frequency at the solar system barycenter using the gtbary utility from the Fermi 
Science Tools package. 

We then searched for coincidences between radio GPs and 7-ray photons. Due to clock 
errors in the backend on four days at 1.2 GHz (55097, 55257, 55290, 55346), only 75131 
MP and 10771 IP GPs were used. However, all of the 5232 MP and 888 IP GPs recorded 
at 330 MHz were used. This resulted in a total of 92022 GPs and 393 7-ray photons, with 
an aver age 7-ray photon ra te of ~16 photons/hour, in agreement with the 15 photons/hour 



seen by lAbdo et al.l (l2010l ). Fermi observes the Crab pulsar for about 11 hours per day, 
resulting in about 34 hours of simultaneous observing time. We also made energy cuts and 
searched for coincidences between radio GPs and 7-ray photons above both 500 MeV and 1 
GeV. 

Based on comparisons with correlations of randomized data, we found no significant 
correlation between MP or IP GPs at both 1.2 GHz and 330 MHz and 7-ray photons out 
to a time lag of ±3 x 10^ spin periods. Since there was not necessarily the same number of 
randomized 7-ray photons during an observation as there were real 7-ray photons, we used 
this large time lag to make sure that there were the same number of real and randomized 
correlations once all of the 7-ray photons were included. The mean and standard deviation 
of the randomized correlations were calculated by randomly assigning 7-ray photon arrival 
times within GTIs and correlating them with GP TOAs. Figure M shows the results of 
this correlation analysis. The random correlations are the mean of 10000 trials, and the 
error bars represent one standard deviation. The largest deviation from the mean is a 2.4cr 
anticorrelation between 7-ray photons and post-fiare 1.2 GHz MP GPs at a time lag of 20000 
spin periods, while the largest correlation is 2.1a for the pre-fiare IP at 1.2 GHz at a time 
lag of 200000 spin periods. The maximum correlation/anticorrelation for each frequency at 



|http : //f ermi . gsf c . nasa . gov/cgi-bin/ssc/LAT/LATDataquery . cgi 
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each energy cut is shown in Table [61 For both the MP and IP for each of the three data sets 
(pre-flare 1.2 GHz, 330 MHz, post-flare 1.2 GHz), there are at most two time lags where the 
correlation exceeds 2a. There are a total of 49 time lags included in the correlation, so two 
time lags are only four percent of the total data, while we statistically expect five percent 
of the data to have correlations beyond 2a. Given the low significance and large time lag of 
these results, they do not provide any compelling case for a physical origin. 

Selecting only 7-ray photons above 500 Mev resulted in a total of 119 photons. The 
maximum correlation for this more restricted set was 2.3cr at a time lag of one spin period 
for pre-flare MP at 1.2 GHz. The maximum anticorrelation was 3.5a at a time lag of 20000 
spin periods for the post-flare MP at 1.2 GHz. No more than one time lag for any data set 
was had a correlation over 2a, so those measurements above 2a are statistically insignificant. 
Statistically, only half of a percent of the data should be above 3(j, and two percent of our 
data have a 3.5(j significance. However, this is likely due to the coarseness of our time lags. 
Given more time lags, we would expect this significance to drop below 3(t. 

Only 65 7-ray photons had energies above 1 GeV. The maximum correlation occured 
in the 1.2 GHz pre-fiare MP with a 3.2cr significance at a time lag of one spin period. The 
maximum anticorrelation was 1.2cr, which occured in the 1.2 GHz pre-fiare IP at a time lag 
of 200 spin periods. As with the 500 MeV energy cut, we only see at most one time lag 
with above 2a, which is insignificant. Similarly, the 3.2a" siginificance is likely due to the 
coarseness of our time lags, and would likely drop given more time lags. 



6. Crab Nebula 7-ray Flare 



The r ecent 7-ray flare fro m the Crab Nebula was detected by the AGILE satellite above 
100 MeV (iTavani et al.ll201ll ). Elevated 7-ray flux was observed from MJDs 5 5457—55461 



No y ariations i n pulse shape were f ound at "Y -rav ([Havs et al 



2010f) . X-ray f lTavani et al. 



2011 ). or radio ( Espinoza et al. 2010) energies. Espinoza et al. ( 201ol ) also found no increase 



in pulsed radio flux, glitches, or changes in DM around the date of the flare. Our closest 
observations before the flare were on MJD 55412 at 330 MHz and MJD 55352 at 1.2 GHz, 
and our first observation after the flare was on MJD 55516 at 1.2 GHz. 

We looked for changes in the average pulse profile, CP shape, power-law index, and 7- 
ray correlation in our pre- and post-flare 1.2 GHz data. We found no significant differences 
in the pulse profile shape and the average CP shape (Figure [TOj) before and on four days 
about two months after the flare. The pre-flare MP power-law indices are in the range 
2.1±0.3 to 3.1±0.2, while the post-flare power-law indices for the MP and IP are 2.56±0.05 
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to 2.93±0.05 and 2.4±0.4 to 2.81±0.03, respectively. We were unable to calculate power- 
law indices for the IP pre-flare due to low IP GP statistics. The largest anticorrelation 
between 7-ray photons and GPs was seen in the post-flare data, but it is still on the order 
of correlations/anticorrelations seen in pre-flare data, both at 1.2 GHz and 330 MHz. 



7. Conclusions 

We compared GB43 and GET GPs and found that 3% of MP GPs and 5% of IP GPs 
at 1.2 GHz were simultaneously detected at 8.9 GHz. The probability of either of these 
events occuring by chance is zero. Also, an additional four IP GPs at 1.2 GHz were within 
one ijs of an IP GP at 8.9 GHz, with a zero percent chance probability. This may suggest 
that, although the folded proflles are much different at the two frequencies, the emission 
mechanism is similar. However, the low percentage of low- and high- frequency IP GPs 
within one fis could mean that the IP emission mechanism is different at higher frequencies, 
especially since the high-frequency IP is shifted by 10 degrees. 

Long observations allowed us to collect the largest sample of GPs to date, which we 
then used to calculate power-law indices for fits to amplitude distributions. These power- 
law indices agree with previously published values. A comparison of these indices to other 
source classes shows that GP emission is not exactly related to emission from other source 
classes and most closely matches the emission from magnetars and RRATs. 

We found no significant correlations between GB43 GPs and Fermi 7-ray photons in the 
energy range 0.1—100 GeV. There are only a few correlations/anticorrelations in the MP and 
IP at both frequencies. They occur at different time lags and are within 2.5a of the mean for 
a correlation with randomized data. This suggests that although increased pair production 
in the magnetosphere may contribute to GP occurence, it is not a dominant factor. More 
likely possibilities for GP generation are increased coherence or changes in beaming. 

Multifrequency correlations have been searched for previously. IShearer et al. J2OO3I) 



found a 3% increase in the brightness of optical pulses at the time of GPs and ICoUins et al. 



(I2OI2I ) found a slight correlation between GPs and enhanced optical pulses, which supports 
our ~2(T correlations/anticorrelations in suggesting that there are small fluctuations in the 
magnetospheric particle density during GPs. The ant icorrelat ions, however, would suggest 
that increased particle density would lead to more radio emission and less 7-ray emission, 
which does not make sense in the context of Lyutikov's theory. Since the signiflcance of 
the anticorrelations are on the order of t he sig niflcance of the correlations, we assume that 



these are insigniflcant. iLundgren et al.l ( 119951 ) found that the 7-ray flux does not vary by 
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more than 2.5 times the average flux during a GP. iBilous et al.l (1201 if ) ruled out a strong 
correlation between GPs at 8.9 GHz and 7-ray photons above 100 MeV. They were able to 
put an upper limit on the 7-ray flux during IP GPs of 8—16 times the average pulsed flux, 
suggesting that there still might be a slight correlation between GPs and 7-ray photons. We 
could not carry out the same experiment because we did not have sufficient 7-ray photons 
to compare the 7-ray profile made with 7-ray photons around GPs with the 7-ray profile 
excluding 7-ray photons around GPs. 



The recent 7-ray fiare from the Crab Nebula occured during the span of our observations, 
so we were able to compare the behavior of the Crab pulsar before and after the fiare. We 
found no significant changes in pulse shape, power-law index, or 7-ray correlation, suggesting 
that there was no change in the properties of the pulsar during the flare. 
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Fig. 1. — Folded pulse profiles from MJD 55406 for Fermi (top), the GET (middle, solid 
line), and the GB43 (bottom, sohd line), as well as histograms of the number of GPs for 
the GET (middle, dashed hue) and GE43 (bottom, dashed hue). The Fermi profile is made 
using photons above 100 MeV over a 24 hour period (~11 hours on source time). The GE43 
observations were taken at a center frequency of 330 MHz over a 220 MHz band for 3.7 hours 
and the GET observations were taken at a center frequency of 350 MHz over a 100 MHz 
band for 10 minutes. The folded radio profiles have been normalized to have peaks of unity. 



This preprint was prepared with the A AS macros v5.2. 
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Fig. 2. — Number of GPs versus pulse phase for all GPs collected with the GB43 at 1.2 GHz 
(top, post- flare; middle, pre- flare) and 330 MHz (bottom). Most GPs at frequencies below 5 
GHz come at the phase (~0) of the MP, while there are still a considerable number at the IP 
phase (~0.4). No GPs are seen at other phases. The FWHMs of the 1.2 GHz distributions 

(226 /is pre-flare, 292 post-fiare) are on the order of the FWHMs of the folded profiles at 
1.2 GHz, while the FWHM of the 330 MHz distribution (704 /is) is much narrower than the 
folded profile and is on the order of the average GP profile at 330 MHz. Even though there 
were only four observing epochs after the flare, there were many more GPs recorded than 
pre-flare due to a large increase in receiver sensitivity. 
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Fig. 3.— A single GP observed both with the GB43 (dashed hne) at 330 MHz and the GET 
(sohd hne) at 350 MHz on MJD 55406. When corrected for the dispersion delay due to 
different observing frequencies, the peaks of the GP from both observations are within one 
fjLS of each other. The GB43 profile has been shifted upwards for clarity. 
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Fig. 4. — Average pulse profiles from all 1.2 GHz (top) and 330 MHz (bottom) observations 
with the GB43, listed in Tabled] The 1.2 GHz profile excludes MJDs 55097, 55257, 55290, 
55299, and 55346, where the folded pulse profiles were dominated by RFI. At 1.2 GHz the 
MP (larger) and IP (smaller) are apparent, an d a low-frequency compo nent can be seen 



one-tenth of a pulse phase ahead of the MP (e.g. iMoffett &: Hankinslll996l ). At 330 MHz the 
MP and IP are much wider due to large DM smearing, which at the top of the band (440 
MHz ] is 0.3 ms and at th e bottom (220 MHz) is 2.4 ms. The scattering time, estimated 
from lKuz'min et al.l ( 120111 ). is 470 fis at the top of the band and 5.8 ms at the bottom, and 
the sampling time is 819.2 /is. This leads to an effective resolution of 990 /is at the top of 
the band and 6.3 ms at the bottom. 
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Fig. 5. — Average GP profiles for all MP (left) and IP (right) GPs (dashed fines) collected 
with the GB43 at 1.2 GHz (top) and 330 MHz (bottom). The 1.2 GHz profiles are made by 
summing the 76707 individual MP and 10871 individual IP GP profiles, and the 330 MHz 
profiles are made by summing the 5232 individual MP and 888 individual IP GP profiles. 
They are shown with the folded MP and IP profiles for all 1.2 GHz and 330 MHz data (solid 
lines). The FWHM of the 1.2 GHz folded MP and IP are ~305 ^s and ~360 ^s, respectively, 
while the FWHM of the average MP and IP GPs are both ~140 fis. The FWHM of the 
330 MHz folded MP and IP are ~3.2 ms and ~1.5 ms, respectively, while the FWHM of the 
average MP and IP GPs are ~1.1 ms. The 330 MHz MP is much wider than the IP due 
to its overlap with the MP precursor. The intensities of the profiles are arbitrary and are 
scaled for clarity. 
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Fig. 6. — A plot of the log of the number of GPs collected with the GB43 vs the log of 
their S/N for the MP at 1.2 GHz on MJD 55099 and a best-fit power- law, which has a slope 
of 2.93±0.07. The deviation from the power-law distribution at high S/N i s seen on all 



days, and these GPs may be examples of supergiant pulses, seen previously by ICordes et al. 



(I2OO4J ). The average supergiant pulse is about 15 times stronger than the average GP. 
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Fig. 7. — Average profiles from tlie GBT at 8.9 GHz (solid line, from iBilous et al.l (l201l[ )) 
and the GB43 at 1.2 GHz (dashed line). The weakening of the MP at high frequencies can 
be seen, as well as the strengthening and shift of the IP. 
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Fig. 8. — Comparison of the strength of the MP of the folded profile to that of the IP for 
each observation. There is an obvious linear trend in both the 330 MHz and 1.2 GHz data, 
with slopes of 0.58±0.08 and 0.310±0.004, respectively. This supports the theory that the 
IP may be a reflection of the MP off of the magnetosphere. 
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Fig. 9. — Number of coincidences between Fermi piiotons and GB43 GP TOAs (solid lines), 
as well as the number of coincidences between randomized Fermi photons and GP TOAs 
(dotted line) with la error bars. The bottom row is the correlation with the 1.2 GHz pre- 
flare data, the middle row is the 330 MHz data, and the top row is the 1.2 GHz post-flare 
data. For each row the MP is the left panel and the IP is the right panel. 
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Fig. 10. — A 7-ray flare occured from MJDs 55457—55461, so we checked for changes in 
the average pulse shape and GP shape. Left: The folded pulse profiles from 1.2 GHz GB43 
observations on MJD 55352 (top) and MJD 55532 (middle), and the difference between them 
(bottom) . The profiles are scaled so that the area under the main pulse for both profiles is 
the same. As can be seen from the bottom plot, the difference is on the order of the rms 
noise, so there is no statistically significant change in pulse shape. Right: The average GP 
profiles from 1.2 GHz GB43 observations on MJD 55352 (top) and MJD 55532 (middle), and 
the difference between them (bottom). This difference is not signiflcant as there is a larger 
variation in average GP shape between prc-flare days. The profiles are also scaled so that 
the area under the main pulse for both profiles is the same. 
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Table 1: Observational parameters including the MJD the observation was taken, the center 
frequency used, the length of the observation, the sampling time, the DM used for dedis- 
persion, and the number of GPs found with S/N > 10. The starred MJDs are days in 
which simultaneous observations with the GBT occured. The dashed line separates pre- and 
post-7-ray flare observations. 
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Table 2: Sampling times, DM smearing times, scattering times, and effective resolution at 
the top and bottom of the band for the GB43 330 MHz and 1.2 GHz observations, as 
well as for the GET 350 MHz observation. The DM smearing times were calculated using 
the delay between t wo consecutive channe ls caused by dispersion, and the scattering times 
were estimated from lKuz'min et al.l ( 120111 ). The sampling time used for the 1.2 GHz GB43 
observation was the most common sampling time from Table [1] 
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Scattering Time 


Effective Resolution 


(MHz) 


(us) 


Top 


Bottom 


Top Bottom 


Top 


Bottom 


330 (GB43) 


819.2 


300 ^s 


2.4 ms 


470 /is 5.8 ms 


990 /is 


6.3 ms 


350 (GET) 


81.92 


361 fis 


855 /is 


670 /is 1.9 ms 


765 /is 


2.0 ms 


1200 (GE43) 


204.8 


23 |US 


180 fis 


4 fis 54 fis 


206 fis 


278 /is 
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Table 3: Differential power-law measurements for the MP and IP for each observing epoch. 
Some MJDs have two entries hsted because there were two separate observations on those 
days. An X means that there were not enough pulses recorded for the MP or IP on that 
day to fit a power-law to the amplitude distribution. The dashed line separates pre- and 
post-7-ray flare observations. 



MJD 


Differential Power- 
law Index fMP) 


Differential Power- 
law Index flPl 




X 


X 


55079 


2.1±0.3 


X 


55088 


X 


X 


55083 


3.1±0.2 


X 


55085 


X 


X 


55086 


2 70 ±0 06 


X 


55090 


2.8±0.1 


X 


55093 


2.6±0.2 


X 


55095 


2.9±0.2 


X 


55096 


X 


X 


55097 


X 


X 


55099 


2.93±0.07 


X 


55102 


2.92±0.08 


X 


55118 


X 


X 


55118 


X 


X 


55139 


2.61±0.08 


X 


55142 


X 


X 


55172 


X 


X 


55178 


2.81±0.08 


X 


55180 


X 


X 


55240 


X 


X 


55240 


2.1±0.1 


X 


55256 


X 


X 
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MJD Differential Power- Differential Power- 
law Index (MP) law Index (IP) 



55257 2.1±0.1 X 

55261 X X 

55264 X X 

55269 X X 

55270 X X 
55284 2.1±0.1 X 
55290 X X 
55299 X X 
55304 X X 

55346 2.4±0.2 X 

55347 2.8±0.1 X 
55347 2.40±0.07 X 
55352 2.21±0.07 X 
55403 2.5±0.2 X 

55405 2.95±0.09 2.4±0.2 

55406 2.94±0.05 3.1±0.2 

55411 2.8±0.1 2.9±0.2 

55412 2.7±0.2 X 
555"l6 2~93±6m 2.4±0.4" 
55532 2.56±0.05 2.7±0.2 
55541 2.8±0.1 2.81±0.03 
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Table 4: Comparison of differential power-law indices for the MP and IP at both 1.2 GHz 
and 330 MHz between this work and previously published values. 



Frequency 
(MHz) 


Differential Power- 
law Index (MP) 


Differential Power- 
law Index (IP) 


Reference 


112 


3.3^^ 


— 


Smirnova & Loe;vinenko f2009) 


146 


3.5 


3.8 


Arevle & Gower fl972) 


200 


2.T 


— 


Bhat et al. f2007) 


330 


2.5-3.0 


2.4-3.1 


This work 


430 


2.3 


_b 


Cordes et al. f2004) 


600 


3.2 


3.0 


PoDov et al. f2009) 


812 


3.3^^ 




Lunderen et al. fl995l 


1200 


2.7-4.2 


2.6 


PoDov & StaDDers f2007) 


1200 


2.1-3.1 


2.4-2.8 


This work 


1300 


2.3^ 




Bhat et al. f2008) 


1400 


2.8 


3.1 


KaruDDUsamv et al. (2010) 


2100 


3.0^^ 




Zhuravlev et al. f 20111 


4850 


2.8^^ 




PoDov et al. r2008) 



^ MP and IP GPs were combined in these analyses 
^ No measurement was taken for the IP 
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Table 5: Differential power-law indices for different source classes for comparison with GP 
power-laws. 



Source 


Power-law 


Reference 


Class 


Index 




Normal Pulsars 


3.85 


Kramer et al. f2002) 


Magnetars 


2.1-7.7 


Servlak et al. f2009) 


RRATs 


3.0 


Miller et al. fin preD.l 



Table 6: Maximum correlation/anticorrelation for each frequency at each energy cut, and the time lags at which they 
occur. Note that the maximum correlation for the 330 MHz IP with the 500 MeV energy cut is negative. 



Energy 


Frequency 


MP/IP 


Max 


Time Lag 


Max Anti- 


Time Lag 


uut 


/ 1\ ,f TT \ 

(Miizj 




Correlation [a) 


(bpm Periods) 


correlation [a) 


(bpm Periods) 




Pre-flare 1200 


MP 


1 J- .T: 


30000 


-0.9 


30 








1 O 1 


zUUUUU 


— i. / 


oUU 


100 MeV 


330 


MP 


-FO.2 


10000 


-2.1 


200 


TT> 

IP 


30 




1 AAA 

1000 


-1-0.5 


— 1.5 




Post-flare 1200 


MP 


+1.5 


40 


-2.4 


20000 




IP 


—2.3 




+1.9 


500 


O A A A A 

30000 




Pre-flare 1200 


MP 


+2.3 


1 


-0.4 


8000 




TP 






+ i.O 


onnnn 
ZUUUU 


— i.o 


( UU 


500 MeV 


330 


MP 


+0.1 


6000 


-1.7 


200 


IP 


300 


-1.7 


2000 




-0.1 




Post-flare 1200 


MP 


+0.8 


200 


-3.5 


20000 




IP 




+1.0 


200 


-2.7 


30000 




Pre-flare 1200 


MP 


+3.2 


1 


-0.1 


60 




IP 






+2.0 


5000 


-1.2 


200 


1 GeV 


330 


MP 


+0.8 


6000 


-0.9 


10 


IP 








+0.9 


6000 


-1.1 


100 




Post-flare 1200 


MP 


+1.9 


60 


-0.4 


1 




10 




IP 


+2.1 


3000 


-0.7 



